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$E 3E. 包 气 带 土壤 水 分 人 渗 在 水 文 循环 中 起 着 至 关 重 要 的 作用 ,深入 了 解 降雨 - 积 水 - 人 渗 过 程 的 相互 关系 对 评 


佑 荒漠 环境 下 土壤 水 分 补给 和 降水 再 分 配 具 有 重要 意义 。 本 研究 以 自然 沙丘 的 丘 间 低 地 为 研究 对 象 ,通过 土壤 剖 
看 入 渗 原 位 观测 试验 。 采 用 Kostiakov , Green-Ampt Philip A 32:78 ffl Hydrus-1D .Hydrus-2D/3D 水 文 模型 模拟 了 上 廿 


壤 水 分 垂直 人 渗 过 程 , 旨 在 寻求 土壤 积 水 人 渗 可 接受 的 方法 。 模 拟 与 实测 结果 比较 表明 :综合 考虑 平方 和 误差 5] 
方 根 误差 等 验证 指标 ,Philip 模型 能 够 预测 砂 土 人 渗 率 .累积 入 渗 量 和 湿润 锋 推进 过 程 ,并且 整个 湿润 区 , Hydrus- 
3D 模拟 土壤 含水 量 的 效果 明显 优 于 Hydrus-2D,RMSE 均值 和 尼 系 数 分 别 达 到 了 0.018 em! * em ”和 0.93。 综 上 所 
示 ,Philip Ej Hydrus-3D 模型 结合 可 有 效 地 估算 砂 士 积 水 人 渗 参 数 和 模拟 水 分 传输 过 程 ,进而 提高 土壤 水 分 人 渗 研 


究 效 率 。 
关键 词 : 包 气 带 ; 积 水 人 渗 ; 数值 模拟 ; 砂 质 土壤 


入 渗 是 指 水 分 进入 土壤 的 过 程 , 而 土壤 水 分 人 
渗 作 为 水 文 循环 中 关键 环节 ,对 连接 地 表 水 资源 和 
降雨 以 及 调控 水 分 分 配 至 关 重 要 '", 许 多 学 者 基于 
不 同 尺度 .不 同 土壤 和 气候 区 域 进 行 了 深入 人 研究。 
在 韭 均 质 土壤 中 提出 了 用 以 表征 入 渗 速 率 、 累 积 入 
渗 量 等 参数 和 人 渗 模 型 ,包括 Kostiakov fll Horton 经 验 
模型 ,Philip 近似 物理 模型 和 Green-Ampt、 基 于 Rich- 
ard 方 程 算法 的 Hydrus-1D/2D/3D 等 物理 模型 。 但 
不 同 土壤 条 件 中 模型 应 用 范围 存在 差异 ,Kostiakov 
方程 简单 ,能 够 很 好 地 拟 合 实测 入 渗 数 据 , 但 其 参数 
无 物理 意义 ,无 法 解析 入 渗 过 程 详细 信息 ,Philip 
是 估算 最 终 人 渗 速 率 较 常 用 的 模型 ,已 在 砂 质 土 
二 "和 寿 土 -壤土 农田 中 得 到 了 验证 ”YY。Green-Ampt 
模型 主要 用 于 均匀 土壤 介质 中 的 积 水 入 渗 ,借助 其 
简便 性 和 物理 基础 ,扩展 模拟 了 一 系列 条 件 下 的 入 
A .降雨 和 径流 过 程 “”。 非 饱和 水 分 运 移 过 程 模 
拟 主要 采用 基于 质量 守恒 定律 和 达 西 定律 的 Rich- 
ards 方 程 来 表达 ,但 Richards 方程 通常 用 非 线性 
描述 ,无 法 求解 析 解 ,导致 计算 过 程 较 复杂 ,特别 在 
复杂 的 初始 和 边界 条 件 中 。 由 美国 农业 部 盐 省 土 
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实验 室 开发 的 Hydrus-1D 软件 可 广泛 用 于 模拟 他 
和 -=- 非 饱和 介质 中 一 维 土 壤 水 分 运 移 ,其 Hydrus- 
2D/3D 软件 应 用 于 模拟 二 维 或 三 维 土壤 水 分 运 移 
和 空间 分 布 , 并 已 获得 大 量 的 验证 和 应 用 1。 

于 旱 区 绿洲 荒漠 区 有 许多 代表 性 的 景观 单元 
类 型 WE WERKA, Sb ERHE 
间 低 地 在 荒漠 生态 系统 水 力 连接 和 水 循环 中 扮演 
着 重要 的 角色 , 事 关 荒漠 生态 系统 的 结构 功能 演 
化 和 绿洲 稳定 。 在 年 均 降 水 量 小 于 200 mm 干旱 区 ， 
降水 具有 降水 稀少 BRR 历时 得 等 特征 ” ,尽管 
降水 量 少 ,但 确 是 干旱 生态 系统 中 唯一 水 分 补给 来 
源 ,尤其 是 储存 在 浅 层 土壤 水 分 调控 地 表 径 流 和 降 
WA RA ia Fe, KP ABE A ERK 
文 过 程 重要 组 成 部 分 ,许多 学 者 已 针对 不 同 尺 度 的 
干旱 荒漠 生态 系统 土壤 包 气 带 水 分 人 渗 方 面 开 展 
了 广泛 研究 ”。 积 水 人 渗 是 由 于 和 渗 强 度 小 于 降 
雨 强度 导致 一 种 压力 人 渗 汪 ,尤其 在 土壤 结 皮 较 厚 
的 干旱 区 ,水 分 停留 时 间 延 长 ,相应 人 渗 量 会 降低 ， 
易 导 致 积 水 入 渗 现 象 发 生 , 而 在 丘 间 低 地 入 渗 速 率 
降低 ,会 导致 下 热 面 土壤 含水 率 提高 1。 
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内 外 有 关 土 壤 入 渗 过 程 研究 方法 较 多 ,如 土 。 壤 水 分 。 


柱 回 填 法 、 单 双环 入 渗 Hood 土壤 入 渗 仪 等 ,其 中 采 
用 单 环 定 水 头 对 水 分 人 渗 进 行 原 位 测量 ,可 用 于 佑 
算 土 壤 水 分 人 渗 特 性 汪汪。 然而 ,由 于 野外 试验 需 
要 花费 大 量 的 时 间 和 费用 , 而 简单 可靠 的 渗透 试 
验 数据 处 理 方 法 显得 尤为 上 必要。 但 由 于 土壤 异 质 
性 ,导致 大 量 入 渗 模 型 在 模拟 土壤 水 分 变化 时 仍然 
有 很 大 局 限 性 。 因 此 ,不 同时 空 条 件 下 ,土壤 水 分 
运 移 过 程 模型 的 识别 和 选择 仍 需 不 断 尝 试 ,基于 


12 原 位 试验 与 设计 

在 选 定 的 试验 区 开展 单 环 入 渗 试 验 , 首 先 去 除 
土壤 表面 可 能 阻碍 入 渗 环 插入 的 残留 物 ,然后 用 落 
锤 将 铁 环 夯实 , 鞠 入 深度 约 为 5 cm( 铁 环 尺 寸 :直径 
20 cm, 高 度 为 30 cm)。 此 外 ,要 注意 保持 铁 环 边缘 
与 土壤 表面 垂直 ,并 尽量 减少 对 土壤 扰动 ,试验 开 
始 前 在 环 内 垫上 滤纸 ,以 防止 水 流 破 坏 土壤 表面 ， 
造成 水 头 不 均匀 (图 2) ,并 沿 着 入 渗 环 下 端 中 心 位 


此 ,本 研究 通过 野外 原 位 入 渗 试验 ,分 析 砂 质 土壤 
中 水 流入 渗 过 程 并 评价 不 同 模 型 在 土壤 水 分 入 渗 
模拟 中 的 性 能 , 旨 在 寻找 降雨 条 件 下 适宜 模型 ,以 
便 合理 预测 沙丘 积 水 入 渗 过 程 , 提 升 对 荒漠 沙丘 局 
地 土壤 包 气 带 水 分 运 移 及 交换 过 程 认 知 。 


1 材料 与 方法 


1.1 研究 区 概况 

研究 在 巴 丹 吉林 沙漠 北部 边缘 的 中 国 科 学 院 
临 泽 内 陆 河流 域 研 究 站 进行 ,该 站 属于 中 国生 态 系 
统 研 究 网 络 (CERN ,地 理 位 置 位 于 100°9'30"E， 
39°24'41"N ,海拔 1381 m。 年 均 气 温 7.5 % ,最 高 气 
温 39 MC, 最 低 气温 -27 % ,年 均 潜在 蒸发 量变 化 为 
1900-2088 mm, 全 年 湿度 变化 为 7.3%~80.9%。 多 
年 平均 降雨 量 为 120 mm, 高 峰 出 现在 7 一 9 月 。 降 
雨 是 研究 区 域 最 敏感 的 气象 因子 ,2006 一 2017 年 该 
区 域 共 接 收 到 降雨 事件 约 500 次 (图 1) ,地 下 水 埋 深 
约 为 4.9 m, 地 下 水 的 饱和 毛管 上 升 不 会 影响 表层 土 
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置 插入 土壤 水 分 观测 探头 PRA ELA Tee T ERE 
依次 为 10 cem、20 cm、40 cm、60 cm,80 cm, Zt fao 
集 顺 连续 收集 土壤 水 分 变化 值 ,时 间 间 隔 1 min 
(图 2)。 
1.3 土壤 样品 收集 与 分 析 

在 研究 区 0~80 cm 内 用 土 外 分 层 采集 原状 土壤 
样品 , 土 样 均 来 自 人 渗 剂 面 中 心 ,然后 烘 干 法 求 土 
BERGE ,浸泡 法 测 饱 和 含水 量 ,激光 粒度 仪 获取 土 
壤 粒 径 组 成 , 定 水 头 法 求 饱和 导 水 率 Ks ,土壤 水 力 
参数 是 水 分 运 移 模拟 的 关键 因素 , 主要 特征 数据 见 
表 1。 
2 模型 描述 
Hydrus-1D/3D 模型 理论 
三 维 水 流 运动 可 以 用 如 下 Richards 方 程 来 描述 : 
= d [emt rw Blew]-s co 


式 中 :0 为 体积 含水 率 (cm + em?) 5h AURA KK 


2.1 


(b) 沙丘 地 貌 


—- 


400 


图 1 试验 区 多 年 降雨 状况 和 干 量 区 沙丘 - 丘 间 低 地 地 貌 


Fig. 1 Location of experimental site, rainfall and image of a mature dune slack in arid region 
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(b) 试验 装置 示意 图 


图 2 单 环 试验 示意 岁 及 土壤 包 气 带 剖面 和 探头 深度 


Fig.2 The profiles of soil vadose zone and schematic diagram of single ring experiment 


表 1 试验 场地 土壤 水 力学 和 物理 特性 


Tab.1 Soil hydraulic and structural properties in experimental site 


土壤 深度 


初始 含水 量 土壤 粒 径 分 布 /9 


土壤 孔 院 


/cm P /cm ecm) 黏 士 (<0.002 mm) ” 粉 填 (0.002 ~ 0.05 mm) 砂 土 (>0.05 mm) 度 /% Madii 
10 1.45 0.04 1.73 10.01 88.26 45.28 wt 
20 1.40 0.01 1.28 74 91.63 46.42 wt 
40 1.42 0.02 1.28 6.14 92.42 46.42 砂 土 
60 141 0.02 0.59 4.53 94.88 46.79 砂 土 
80 1.43 0.03 0.79 5.66 93.55 46.04 砂 土 


Cem) ;1 为 时 间 (h);K(h) 为 非 饱 和 导 水 率 函 数 (cm-: 
d);z 为 向 上 正 的 空间 坐标 (m) ;5 为 汇源 项 ,通常 表 
AN AR BAK (S) ,本 试验 在 裸露 沙丘 进行 ,不 考虑 
根系 吸水 。 其 中 简化 的 一 维 水 流 运动 方程 如 下 : 


aog o 


2.1.1 土壤 水 力 性 质 用 Van Genuchten- Mualem 
(VG) 模 型 撕 述 土壤 水 分 特征 曲线 ,其 土壤 水 力 函 数 
如 下 : 


| o 
K(e) - K,S! í - [ - si] (4) 


式 中 :0. 和 6. 分 别 为 土壤 残余 含水 量 和 土壤 饱和 含水 
Ht Cem? * em?) ; S.A AEE (Dot 290) ;kK 为 饱和 导 水 
率 (cmd');n 为 孔径 分 布 指数 (无 量 纲 ) 5 m-1- Vn. 
为 土壤 水 分 特征 曲线 参数 ;gq 为 土壤 基质 势 ;1 为 筷 


际 连 通 性 参数 (无 量 纲 )。 基 于 VG 模型 土壤 水 力 参 
数 见 表 2。 
2.1.2 初始 和 边界 条 件 

(1) 基于 Hydrus-2D 二 维 计算 模拟 

模拟 区 域 人 渗 面 以 上 DA BC 为 不 透水 边界 ;地 
A DE FC 为 大 气 边界 ,下 底面 AB 为 自由 排水 边界 ; 
FEF 在 入 渗 开始 后 很 快 达到 饱和 ,为 定 水 头 边界 ,其 
边界 和 初始 条 件 如 下 : 

CD 初始 条 件 : 

0(x,y,0) -0,,0 cmxx«200 cm; 0 cms y< 100 cm; 
t=0 
式 中 : 0, WI 7K et Cem? em?) 

@ 边界 条 件 : 

DE、FC 边界 :- Kma + K(h) = E(), 0 em<x< 


200 cm; y=100 cm; t > 0 
DA 、BC 边界:- Kma =0, 0 cmsxs 200 cm; 
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表 2 试验 地 不 同 剖 面 土壤 水 力 人 参数 组 成 


Tab.2 Soil parameters for different soil layers in the experiment site 


土壤 深度 /em (UA KE em" em?) RAR AACE m em?) — 3E fem 孔径 分 布 指数 MATKE em: d) 
10 0.373 0.062 0.012 2.292 26.85 
20 0.376 0.071 0.014 2.600 38.41 
40 0.363 0.059 0.013 2.301 40.28 
60 0.358 0.061 0.01 2.475 40.54 
80 0.375 0.082 0.016 2.145 16.76 
0 ems ys 100 cm 0 , x=0 cm 2K 200 cm; -KO =0 em, z=0 cm 或 100 
JR. 28. ROS a ki - 
FI :0 (0, Y t)= 0., K(h) ay + K(h) H() t> cm 
0; » GHJ 边界 : -K(n)% +K (h) +K(h) =0, (50; 
边界 : -KO h $ K(h) =0, 1>0; 0 em< x< 200 0 cm< x< 200 cm; 0 em<z< 100 cm。 
em; y=0 二 维 与 三 维 模拟 均 采用 1 em 压力 水 头 作为 单 


RP: BO 为 人 渗 速 率 ; HO) ATA AK 

(2) 基于 Hydrus-3D 三 维 计 算 模拟 

三 维 模拟 区 域 人 渗 面 OPNM 为 大 气 边 界 ;侧面 
OPGH , GPMJ , MNIJ 和 HION 为 不 透水 边界 ,下 底面 
GHJIJ 为 自由 排水 边界 ;以 直径 R=20 cm 入 渗 圆 面 开 
始 后 很 快 达到 饱和 ,为 定 水 头 边 界 ,其 边界 和 初始 
条 件 如 下 : 

CD 初始 条 件 : 

A(x, Y, Z, 0) =00, 0 cm<x<200 cm; 0 cm< y< 100 
cm; 0 cms z< 100 cm; t=0 

(2 边界 条 件 : 

OPNM 边界 : -Kh ðh 20. 0 em<x<200 cm; 0 


cm&z& 100 cm; y=100 cm; M 
OPGH , GPMJ, MNIJ 和 HION 边界 : -Kht 


by K DD VUUVVVNNNNS 


o i32 


ALIQ qb 
AY KOR UU TERK 


A " E Ek E 


不 透水 边界 
不 透水 边界 


环 和 人 渗 下 的 水 流 边界 条 件 ( 图 3) ,并 将 模拟 区 域 离 
散 成 径 向 网 格 间 距 为 2 em 垂直 网 格 间距 为 1 em 的 
有 限 元 ,总 模拟 时 长 12h。 

22 入 渗 模 型 理论 

2.2.1 Green-Ampt Green-Ampt 模 型 最 初 是 利用 达 


西 定律 对 定 水 头条 件 下 土 柱 积 水 入 渗 进 行 分 析 , 入 
渗 方 程 如 下 : 
H,+S, 
Rd (5) 


ZH S. 
= eut -x 


s Z; 
Z= (6) 
0,65 
x 5j 4 (7) 


SUA SIT BARA ft (em) Lg AS XE (em min); 
KK, 为 饱和 导 水 率 ;5 为 湿润 锋面 吸力 ;如 为 概 化 的 湿 


p ”大 气 边界 M 
R-20 cm 


Ti sg zK3k-1 em 


不 透水 边界 


图 3 Hydrus-2D/3D 模 拟 区 域 及 边界 条 件 
Fig.3 The Hydrus-2D/3D simulation domain and boundary conditions 
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润 锋 深度 (cm) ;6 为 初始 土壤 含水 量 (cm ecm)。 拟 合 度 越 好 。 
分 层 均 质 土壤 入 渗 时 ,通过 整合 Green-Ampt 模 
型 推导 出 累积 人 渗 量 和 湿润 锋 随 时 间 变 化 的 隐 函 3 结果 与 分 析 


数 ,分 别 表示 如 下 : 
/= Yo. ~ 6,;) * (z, ~ yD 851) (8) 


StH: DALEIBE (em); Pius JN aa RS 
始 状态 .饱和 状态 EJ AU 
— O vat = Ori 

N K. 


Z,*S,*H (9) 
e PE Te] 
式 中 ;为 湿润 锋 到 达 第 N 层 界面 所 需 时 间 (min)。 
2.2.2 Philip Philip 模型 最 初 由 Philip 提出 ,累积 
渗 量 和 入 渗 率 表示 为 : 
i(1) 20.5807 +A (10) 
I(t) 2 Si^ At (11) 
式 中 ;4 是 关于 土壤 性 质 和 水 分 导 水 系数 函数 (em- 
min!) ;5 为 吸水 率 ,是 土壤 基质 势 的 数 (cm:min”)。 
ZA0=B - S.A (12) 
式 中 :A9 为 饱和 含水 量 和 初始 含水 量 的 差 值 。 
2.23 Kostiakov-Lewis 修正 后 长 周期 的 Kostiakov 模 
Ak. 


i(t)=Br +i, (13) 


I()= 1E 6 eig (14) 


式 中 :B.C 分 别 为 方程 参数 (C>0, 0<B<1) i) 分 
别 为 人 渗 率 和 稳定 入 渗 率 。 
2.3 模拟 效果 评价 

模拟 值 与 观测 值 比较 用 以 评价 模型 性 能 ,本 研 
究 以 决定 系数 (RR) 和 均 方 根 误差 (RMSE ) 作 为 最 终 
评价 指标 。 


R=1- pa (Xoo: = E (15) 
E NE E 
RMSE- aos -X noli) (16) 


式 中 ;Xi 为 观测 值 ;Xi 为 模拟 值 ; Xaa 为 观测 样 
本 均值 ,obs, i 代表 第 i(1,2,…,n) 个 观测 值 ;n AW 
测 数据 总 个 数 。RMSE 指标 值 越 小 ,表明 模拟 误差 
越 小 ,而 尼 的 指标 值 越 接近 1, 表明 模拟 精度 越 高 ， 


3.1 实测 数据 

结果 表明 ,试验 开始 后 入 渗 速 率 持续 下 降 ,300 
min 左右 达到 稳定 人 渗 速 率 0.57 cm.min (图 4a) , 
累积 入 渗 量 随时 间 增 加 而 递增 ,观测 时 上 段 内 累积 
渗 量 达到 398 cm( 图 4b), 入 渗 开 始 后 湿润 锋 同时 沿 
竖 直 和 水 平 两 个 方向 快速 推进 ,200 min 后 推进 速度 
趋 缓 ,入 渗 结 束 后 ,水 平方 向 湿润 锋 距 离 较 垂直 方 
向 湿润 锋 距 离 增加 了 45% (图 4c)。 可 以 发 现 震 函 
BAT Vete m HUE SU HE (R°=0.85) .. ZA TIU , AF PAI 
数 只 是 估算 土壤 水 分 人 渗 的 一 个 简单 的 经 验 模型 ， 
无 具体 物理 参数 意义 。 此 外 ,结果 表明 ,在 10~80 
cm 深度 内 土壤 各 个 剖面 含水 量 差异 明显 (图 4d), 达 
到 稳 态 土壤 含水 量 随 土 层 深度 增加 而 降低 ,但 20 
cm 剖面 的 最 大 含水 率 明 显 低 于 其 他 剖面 , 且 均 低 于 
饱和 含水 量 ,这 可 能 与 容重 REA AE RETE 
差异 有 关 。 
3.2 观测 与 模拟 结果 比较 
32.1 ABRE 4 种 水 文 模型 的 入 渗 率 模拟 值 与 
观测 值 对 比如 图 5 所 示 。 结 果 表 明 , 入 渗 率 的 让 系 
数 均值 在 0.75~0.95, Philip 的 尼 系 数 最 高 (0.95) ， 
Kostiakov 的 模拟 结果 略 大 于 观测 值 ,相反 Green- 
Ampt 的 模拟 结果 略 低 于 观测 值 , 尤 其 在 人 渗 后 期 。 
Hydrus- 1D 模拟 结果 与 实测 匹配 较 低 ,决定 系数 尼 
较 小 (0.79)。 总 之 , 芒 漠 砂 质 十 壤 积 水 人 活 条 件 下 ， 
Philip 较 其 他 入 渗 模 型 能 更 好 地 描述 土壤 水 分 人 渗 
3.2.2 READE 模拟 与 观测 累积 入 渗 量 关系 如 
图 6 所 示 ,结果 表明 ,4 种 模型 累积 入 渗 量 模拟 值 均 
与 实测 数据 有 较 好 的 相关 性 。 而 Philip 模型 尼 最 
高 ,RMSE 系数 最 低 ,Kostiakov 模型 RMSE 次 之 ( 表 
3) ,但 Hydrus-1D 模型 下 RMSE 值 系数 相对 较 高 , 表 
明 高 估 了 累积 入 渗 量 ,尤其 在 入 渗 初 期 ,但 4 种 入 渗 
模型 尼 系 数 均 达 到 了 0.98 以 上 ,因此 ,4 种 人 渗 模 型 


均 是 预测 砂 士 累积 人 渗 量 较为 有 效 的 模型 ,尤其 是 


Philip 模型 。 

3.2.3 湿润 锋 距 离 ”模拟 与 观测 湿润 锋 距 离 之 间 的 
关系 如 图 7 所 示 ,结果 表明 ,Green-Ampt 和 Hydrus- 
ID 模型 的 湿润 锋 距离 模拟 值 与 观测 值 差异 较 大 , 尤 
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0.75 F (a) 入 渗 速 率 随 时 间 变化 (b) 累积 入 渗 量 随时 间 变 化 
Bp 
t 070L 
E e 观测 结果 a 
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图 4 观测 时 段 的 入 渗 速 率 、 累 积 入 渗 量 .湿润 锋 距 离 和 土壤 含水 量变 化 
Fig. 4 Dynamic of observed infiltration rate, cumulative infiltration, depth of wetting front and 
soil water content during the study period 
去 0.80: F (a) Kostiakov 模 型 zx 0:80 (b) Green-Ampt 模 型 
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1 1 1 
05 06 07 08 09 10 
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注 : 直 线 表 示 数 据 之 间 可 能 存在 1:1 关 系 。 下 同 


图 5 入 渗 速 率 模拟 值 与 实测 值 比较 


Fig. 5 Comparison of simulated infiltration rate from models with observed result 
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图 6 累积 入 渗 量 模拟 值 与 实测 值 比较 


Fig.6 Comparison of simulated cumulative infiltration from models with observed result 
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图 7 湿润 锋 距离 模拟 值 与 观测 值 比较 
Fig.7 Comparison of simulated wetting front depth from models with observed result 

其 是 Hydrus-1D 明显 高 舍 了 湿润 锋 推进 距离 ,其 均 。 尺 均 值 结果 可 知 ,Philip 模 型 对 砂 土 湿润 锋 推进 距离 
(E R/F Philip 模型 。 一 种 可 能 的 解释 是 Hydrus- ”的 预测 效果 较 好 。 

ID 忽略 了 土壤 水 侧 渗 和 气流, 以往 研究 表明 ,Hy- 32.4 土壤 含水 量 与 Hydrus-2D 模拟 Æ F Hy- 
drus-1D 高 估 了 湿润 区 的 蓄 水 能 力 , 但 综合 RMSE 和 ”drus-2D 模 拟 土壤 含水 量 如 图 8 所 示 。 结 果 表 明 , 土 


202201.00111v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


130 干旱 区 研究 39 卷 
表 3 模拟 结果 拟 合 优 度 参数 
Tab.3 Goodness-of-fit parameters for simulation results with models 
E Kostiakov Green-Ampt Philip Hydrus-1D 
试验 测定 项 目 2 
RMSE/(cm* cm?) R RMSE/(cm*+ cm”) R RMSE/(cm?* cm?) R RMSE/(cm*+ cm?) R 
AGB EXE (cm * min?) 0.003 0.85 0.002 0.84 0.0003 0.95 0.004 075 
累积 人 渗 量 /cm 0.08 0.98 0.35 0.97 0.04 0.99 1.38 0.98 
湿润 锋 距 离 /cm - - 1.85 0.81 0.24 0.94 - 0.27 
IE: RMSE 表示 均 方 根 误差 , 尼 表 示 决 定 系 数 。 下 同 。 
_ (a) 10 cm (c) 40 cm 
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图 8 基于 Hydrus-2D 模 拟 土壤 含水 量 与 实测 值 比较 


Fig. 8 Comparison of simulated soil water content by Hydrus-2D model with observed result 


壤 水 分 模拟 值 与 观测 值 拟 合 度 较 低 ,RR 均值 仅 为 
0.82 ,尤其 在 表层 0~10 em, RAH 0.65 ,综合 0~80 cm 
的 RMSE 和 尼 均 值 结 果 可 知 ,Hydrus-2D 模型 并 不 能 


很 好 地 预测 积 水 人 渗 条 件 下 土壤 水 分 变化 过 程 。 
32.5 土壤 含水 量 与 Hydrus-3D 模拟 基于 Hy- 
drus-3D 模拟 土壤 含水 量 如 图 9 所 示 ,结果 表明 , 土 
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图 9 基于 Hydrus-3D 模 拟 土壤 含水 量 与 实测 值 比 较 


Fig.9 Comparison of simulated soil water content by Hydrus-3D model with observed result 


壤 含 水 量 模拟 值 与 观测 值 基本 一 致 ,0~80 cm 土壤 
剖面 的 尺 均值 为 0.93,RMSE 均值 为 0.02 em? .cm 
( 表 4)。 尽 管 较 Hydrus-2D 4H EE , Hydrus-3D 对 入 
渗 后 期 土壤 含水 量 预 测 较 低 ,但 是 综合 考虑 尼 和 
RMSE 值 ,Hydrus-3D 是 描述 和 模拟 砂 土 积 水 条 件 土 
壤 仿 水量 较 为 理想 的 选择 。 


4 讨论 


41 土壤 水 分 入 渗 特 性 
单 环 积 水 入 渗 仪 是 测量 土壤 水 分 入 渗 速 率 的 
最 常用 方法 之 一 ,已 被 广泛 用 于 评估 和 监测 降雨 人 


渗 进 入 土壤 过 程 。 研 究 表明 , 随 着 时 间 的 延长 ,入 
渗 速 率 随 着 人 渗 量 的 增加 而 降低 ,最 终 达到 稳定 人 
渗 率 ””。 本 研究 发 现 , 砂 土 的 稳定 入 渗 率 为 0.55 
cm'min ,这 与 前 人 研究 砂 士 和 人 渗 速率 范围 为 0.2~ 
0.77 em: min Zt RAN $7 ,结果 差异 可 能 由 于 土 
二 质地 初始 土壤 含水 量 和 地 表 覆 盖 条 件 导致。 土 
壤 水 分 运 移 过 程 中 , 非 饱和 区 水 分 湿润 锋 推进 动力 
学 特性 研究 具有 重要 的 意义 ,人 研究 指出 , 单 环 入 渗 
受 环 下 水 分 横向 运动 影响 ,其 基质 势 通 量 在 湿润 锋 
处 占 主导 地 位 ”, 本 研究 发 现 湿 润 锋 距离 水 平方 向 
较 垂直 方向 大 。 此 外 ,在 入 渗 发 生 的 前 100 min Pj, 
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#24 H ydrus-2D f Hydrus-3D 模拟 结果 拟 合 度 参数 
Tab.4 Goodness-of-fit parameters for simulation results 
between Hydrus-2D and Hydrus-3D 


Hydrus-2D Hydrus-3D 
土壤 深度 /em RMSE RMSE 
fcm? * cm?) /(cm em?) 
10 0.02 0.65 0.013 0.94 
20 0.024 0.83 0.021 0.84 
40 0.017 0.91 0.026 0.96 
60 0.032 0.85 0.021 0.93 
80 0.015 0.94 0.009 0.98 


湿润 锋 推 进 速度 随时 间 增 加 而 减 小 ,之 后 保持 恒定 
不 变 ( 图 10) ,一 种 可 能 的 解释 是 砂 士 引发 的 漏斗 流 
在 水 平方 向 上 受 基 质 势 驱动 ,而 在 垂直 方向 上 受 重 
力 势 驱动 。 


4.0 

3.5 。 垂直 方向 观测 值 
o 水 平方 向 观测 值 

3.0 fp 一 拟 合 曲线 


湿润 锋 推 进 速 率 /(cmminD 


1 1 
100 200 300 400 500 600 


700 800 
累积 入 渗 时 间 /min 
图 10 垂直 和 水 平 两 个 方向 湿润 锋 推进 速度 


Fig. 10 Comparison of observed the velocity of the wetting 


front with two directions 


先前 研究 表明 , Philip 模型 可 用 于 砂 质 土壤 人 
渗 速 率 与 人 渗 时 间 关 系 预 测 ” ,Zolfaghari 等 也 做 了 
ASW FRI? ,然而 ,Duan 等 人 和 Ogbe 等 通过 人 研究 
人 渗 模 型 对 砂 土 水 分 入 渗 适 应 性 发 现 ,Hortons 模 型 
较 适 宜 累 积 入 渗 量 预测 , Wang 等 ”通过 模拟 与 实测 
值 比较 ,证 实 Hydrus-1D 能 够 预测 层 状 土 柱 砂 土 积 
水 条 件 下 的 入 渗 速 率 。 然 而 ,由 于 模型 参数 和 性 能 
随 土壤 条 件 和 时 间 的 不 同 而 存在 差异 ,需要 大 量 试 
验 以 确定 入 渗 模 型 适用 范围 和 条 件 。 
4.2 土壤 剖面 水 分 状况 

土壤 水 分 分 布 状 况 是 影响 干旱 区 荒漠 植物 根 
系 吸水 有 效 性 的 重要 因素 ,本 人 研究 结果 表明 ,土壤 
含水 量 表层 高 于 深层 , 积 水 人 渗 导致 土壤 水 分 在 整 
个 湿润 区 呈 不 均匀 非 饱和 状态 分 布 , 这 可 能 与 士 层 


结构 对 土壤 的 持 水 性 能 影响 有 关 ,尽管 整个 剖面 土 
壤 质 地 为 砂 土 ,但 由 于 各 层 之 间 粒 径 组 成 比例 有 差 
异 ,导致 持 水 性 咯 有 不 同 ,出现 土壤 含水 量 空间 变 
SR HER , Hydrus-3D ££ Hydrus-2D 对 土壤 含水 
量 数值 模拟 结果 效果 最 佳 ,其 优点 是 考虑 了 水 分 在 
土壤 中 的 三 维 变 饱 和 多 和 孔 介 质 运动 。 


5 结论 


以 自然 沙丘 中 的 丘 间 低 地 为 研究 对 象 ,通过 单 
环 人 渗 试 验 , 以 确定 土壤 水 力 参 数 ,并 分 析 了 人 砂 土 
在 固定 水 头 渗流 区 的 土壤 水 分 变化 对 入 渗 响 应 ,以 
验证 几 种 入 渗 模型 在 积 水 入 渗 条 件 下 的 适应 性 状 
况 , 取 得 以 下 主要 结论 : 

(1) 入 渗 初 期 人 渗 速率 随时 间 急 剧 下 降 , 之 后 
趋 于 稳定 ,整个 人 渗 过 程 中 ,水平 湿 润 锋 的 推进 速 
度 较 垂 直方 向 快 ,入 闯 速 率 、 累 积 人 渗 量 和 湿润 锋 
深度 与 时 间 呈 震 函 数 关系 。 

(2) 受 边界 和 初始 条 件 以 及 土壤 水 力 参数 设置 
等 因素 影响 ,Hydrus-1D K OCBE TE EKIA 8 
模拟 中 表现 欠 佳 ,但 当前 试验 条 件 下 ,通过 性 能 评 
价 指标 系数 比较 ,特定 土壤 环境 中 Philip 较 Green- 
Ampt, Kostiakov- Lewis 和 Hydrus- 1 D 模型 能 够 合理 
描述 人 砂 土 信 渗 参 数 , 反 映 较真 实 的 土壤 水 动力 
过 程 。 

(3) 单 环 入 渗 试验 需要 同时 考虑 环 内 一 维和 环 
外 三 维 水 流 过 程 ,而 Hydrus-3D 模 型 基本 能 反映 砂 
土 不 同 训 面 三 维 土壤 水 分 变化 。 总 之 ,结合 Philip 
和 Hydrus-3D 水 文 模型 确定 人 砂 土 入 渗 是 可 行 的 ,但 
此 结果 只 是 基于 人 砂 土 , 能 否 应 用 其 他 类 型 土壤 水 分 
入 渗 过 程 的 评估 有 符 进 一 步 研究 。 
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A comparative study of ponded infiltration in a desert sandy soil based 
on multi-hydrological models 
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Abstract: Water infiltration of the vadose zone plays a key role in the water cycle. A better understanding of the 
relationships among rainfall, ponded water, and soil infiltration processes is critical to estimate groundwater 
replenishment and redistribution in desert environments. However, research on the moisture distribution in the 
vadose zone of sandy soil is still limited. In particular, few studies have examined the water transport processes in 
the desert vadose zone. Therefore, we aimed to determine the soil moisture distribution and variation in 
homogenous sandy soil under ponded infiltration and to validate models using data from these observations. This 
study was conducted in a nature dune slack through dynamic in situ observations using profile infiltration 
experiment and continuous soil moisture measurements at 1 min intervals. A single-ring infiltrometer (diameter, 
20 cm; insertion depth, 5 cm) was used to keep a constant water head. Kostiakov- Lewis, Green- Ampt, Philip, 
Hydrus-1D/2D/3D numerical software based on Richard' s equation were used to describe the vertical infiltration 
and soil water movement processes and to identify suitable models for ponded infiltration processes in sandy soil. 
The validation indices included the sum of squared error and root mean squ are error. A comparison between the 
simulated and measured results indicated that the Philip model could predict the infiltration rate, cumulative 
infiltration, and wetting front advancement correctly, with an RMSE value of 0.0003, 0.04, and 0.24 cm min ' 
and R^ of 0.95, 0.99, and 0.94, respectively. These values were clearly less than those of the Kostiakov- Lewis, 
Green- Ampt, and Hydrus- 1D models, although the other models also showed good overall agreement with the 
field measured cumulative infiltration. The Hydrus- 3D model yielded a better fit to the simulated soil moisture 
than the Hydrus-2D in the wetted zone based on the coefficient of determination, average RMSE (0.018 cm cm ?) 
and A^ (0.93). Altogether, it appears that the combination of the Philip and Hydrus- 3D models is a highly 
effective approach to estimate ponded infiltration parameters and simulate the water transport process in variably 
saturated sandy soil. The three-dimensional soil water diffusivity must be considered to predict water infiltration 
in sandy soils. Our results also demonstrate that the integration of numerical simulation and field measurements 
can improve the research efficiency of soil water infiltration in dry environments. In practice, the accurate 
prediction of soil water infiltration can be done by the selection of the proper models based on the soil properties 
of the particular sites. 


Keywords: vadose zone; ponded infiltration; numerical simulation; sandy soil 


